Constraints on holographic dark energy from X-ray gas mass fraction of galaxy 

clusters 
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We use the Chandra measurements of the X-ray gas mass fraction of 26 rich clusters released 
by Allen et al. to perform constraints on the holographic dark energy model. The constraints are 
consistent with those from other cosmological tests, especially with the results of a joint analysis 
of supernovae, cosmic microwave background, and large scale structure data. From this test, the 
holographic dark energy also tends to behave as a quintom-type dark energy. 



Recent observations of Type la supernovae (SNe la) 
indicate that the expansion of the Universe is acceler- 
ating at the present time. These results, when combined 
with the observations of cosmic microwave background 
(CMB) and large scale structure (LSS) 0, strongly 
suggest that the Universe is spatially flat and dominated 
by an exotic component with large negative pressure, re- 
ferred to as dark energy |4| . The first year result of the 
Wilkinson Microwave Anisotropy Probe (WMAP) shows 
that dark energy occupies about 73% of the energy of 
our Universe, and dark matter about 23%. The usual 
baryon matter which can be described by our known par- 
ticle theory occupies only about 4% of the total energy 
of the Universe. Although we can affirm that the ulti- 
mate fate of the Universe is determined by the feature 
of dark energy, the nature of dark energy as well as its 
cosmological origin remain enigmatic at present. The 
most obvious theoretical candidate of dark energy is the 
cosmological constant A which has the equation of state 
w = —1. An alternative proposal is the dynamical dark 
energy (quintessence) [B| which suggests that the energy 
form with negative pressure is provided by a scalar field 
evolving down a proper potential. The feature of this 
class of models is that the equation of state of dark en- 
ergy w evolves dynamically during the expansion of the 
Universe. However, as is well known, there are two diffi- 
culties arise from all these scenarios, namely the two dark 
energy (or cosmological constant) problems — the fine- 
tuning problem and the "cosmic coincidence" problem. 
The fine-tuning problem asks why the dark energy den- 
sity today is so small compared to typical particle scales. 
The dark energy density is of order 10 _47 GeV 4 , which 
appears to require the introduction of a new mass scale 
14 or so orders of magnitude smaller than the electroweak 
scale. The second difficulty, the cosmic coincidence prob- 
lem, states "Since the energy densities of dark energy and 
dark matter scale so differently during the expansion of 
the Universe, why are they nearly equal today"? To get 
this coincidence, it appears that their ratio must be set to 
a specific, infinitesimal value in the very early Universe. 

Recently, considerable interest has been stimulated in 
explaining the observed dark energy by the holographic 
dark energy model. For an effective field theory in a 



box of size L, with UV cut-off A c the entropy S scales 
extensively, S ~ L 3 A 3 . However, the peculiar thermo- 
dynamics of black hole jfj has led Bekenstein to postu- 
late that the maximum entropy in a box of volume L 3 
behaves nonextensively, growing only as the area of the 
box, i.e. there is a so-called Bekenstein entropy bound, 
S < Sbh = nAl 2 L 2 . This nonextensive scaling sug- 
gests that quantum field theory breaks down in large 
volume. To reconcile this breakdown with the success 
of local quantum field theory in describing observed par- 
ticle phenomenology, Cohen et al. proposed a more 
restrictive bound - the energy bound. They pointed out 
that in quantum field theory a short distance (UV) cut- 
off is related to a long distance (1R) cut-off due to the 
limit set by forming a black hole. In other words, if the 
quantum zero-point energy density px is relevant to a UV 
cut-off, the total energy of the whole system with size L 
should not exceed the mass of a black hole of the same 
size, thus we have L 3 px < LM 2 . This means that the 
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maximum entropy is in order of Sg H . When we take the 
whole Universe into account, the vacuum energy related 
to this holographic principle [8j is viewed as dark energy, 
usually dubbed holographic dark energy. The largest IR 
cut-off L is chosen by saturating the inequality so that 
we get the holographic dark energy density 



px = 3c 2 M 2 L~ 2 



(1) 



where c is a numerical constant, and M p = 1/^/8ttG is 
the reduced Planck mass. If we take L as the size of the 
current Universe, for instance the Hubble scale if -1 , then 
the dark energy density will be close to the observed data. 
However, Hsu p pointed out that th is y ields a wrong 
equation of state for dark energy. Li [T(| subsequently 
proposed that the IR cut-off L should be taken as the 
size of the future event horizon 
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Then the problem can be solved nicely and the holo- 
graphic dark energy model can thus be constructed suc- 
cessfully. The holographic dark energy scenario may pro- 
vide simultaneously natural solutions to both dark energy 
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problems as demonstrated in Ref-EU- For related work 

see ElElElElElEi 

Consider now a spatially flat FRW (Friedmann- 
Robertson- Walker) Universe with matter component p m 
(including both baryon matter and cold dark matter) and 
holographic dark energy component px, the Friedmann 
equation reads 



3M 2 H 2 



Pm + PX 



or equivalently, 
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Note that we always assume spatial flatness throughout 
this paper as motivated by inflation. Combining the def- 
inition of the holographic dark energy JJ) and the defini- 
tion of the future event horizon ©, we derive 
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We notice that the Friedmann equation (0J implies 
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Substituting © into (0 , one obtains the following equa- 
tion 



y/l - n x dx' = ce x/2 
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where x = In a. Then taking derivative with respect to 
x in both sides of the above relation, we get easily the 
dynamics satisfied by the dark energy, i.e. the differential 
equation about the fractional density of dark energy, 
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(8) 



This equation describes behavior of the holographic dark 
energy completely, and it can be solved exactly 0, 0] . 
From the energy conservation equation of the dark en- 
ergy, the equation of state of the dark energy can be 
given E3 
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Note that the formula px — Pm a 3 an d the differ- 
ential equation of fix © ar e used in the second equal 
sign. It can be seen clearly that the equation of state 
of the holographic dark energy evolves dynamically and 
satisfies -(l+2/c)/3 < w < -1/3 due to < fl x < 1. In 
this sense, this model should be attributed to the class 
of dynamical dark energy models even though without 
quintessence scalar field. The parameter c plays a signifi- 
cant role in this model. If one takes c = 1, the behavior of 
the holographic dark energy will be more and more like a 



cosmological constant with the expansion of the Universe, 
and the ultimate fate of the Universe will be enterin g th e 
de Sitter phase in the far future. As is shown in Ref . [l0| . 
if one puts the parameter £l x — 0.73 into @, then a 
definite prediction of this model, wo = —0.903, will be 
given. On the other hand, if c < 1, the holographic dark 
energy will behave like a quintom-type dark energy pro- 
posed recently in Ref.|17|. the amazing feature of which 
is that the equation of state of dark energy component 
w crosses the phantom divide, —1, i.e. it is larger than 
— 1 in the past while less than —1 near today. The recent 
fits to current SNe la data with parametrization of the 
equation of state of dark energy find that the quintom- 
type dark energy is mildly favored 0, . Usually the 
quintom dark energy model is realized in terms of double 
scalar fields, one is a normal scalar field and the other is 
a phantom-type scalar field 20] (for quintom model see 
e.g. |21|). However, the holographic dark energy in the 
case c < 1 provides us with a more natural realization 
for the quintom picture. If c > 1, the equation of state 
of dark energy will be always larger than —1 such that 
the Universe avoids entering the de Sitter phase and the 
Big Rip phase. Hence, we see explicitly, the determina- 
tion of the value of c is a key point to the feature of the 
holographic dark energy as well as the ultimate fate of 
the Universe. 

The holographic dark energy model has been tested 
and constrained by various astronomical observations 
[l2l IT3L Il6l | . In a recent work Ell , it has been explicitly 
shown that regarding the latest supernova data as well 
as the CMB and LSS data, the holographic dark energy 
behaves like a quintom-type dark energy. This indicates 
that the numerical parameter c in the model is less than 
1. The best fit results provided by ll| are: c = 0.81, 
£!° n = 0.28, and h = 0.65, which lead to the present 
equation of state of dark energy wq — —1.03 and the 
deceleration/acceleration transition redshift zt = 0.63. 
It is necessary to test dark energy models and constrain 
their parameters using as many techniques as possible. 
Different tests might provide different constraints on the 
parameters of the model, and a comparison of results de- 
termined from different methods allows us to make con- 
sistency checks. In this Letter, we use the X-ray gas 
mass fraction of rich clusters, as a function of redshift, 
to constrain the holographic dark energy model, and to 
compare the results with the previous analysis. 

The matter content of the largest clusters of galaxies 
is thought to provide an almost fair sample of the matter 
content of the Universe. A comparison of the gas mass 
fraction of galaxy clusters, / gas = M gas /M tot , inferred 
from X-ray observations, with determined by nucle- 
osynthesis can be used to constrain the density parameter 
of the Universe ft^ directly [13 . Sasaki [H and Pen H] 
were the first to describe how the / gas data of clusters of 
galaxies at different redshifts could also, in principle, be 
used to constrain the geometry and, therefore, dark en- 
ergy relevant parameters of the Universe. The geometri- 
cal constraint arises from the fact that the measured / gas 
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values for each galaxy cluster depend on the assumed an- 



gular diameter distances to the clusters as / g 
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The measure d f<r a, s values should be invariant with red- 
shift [23L I24I l25| when the reference cosmology used in 
making the measurements matches the true, underlying 
cosmology. The first successful application of such a test 
to constrain cosmological parameters was carried out by 
Allen et al. (2(|; see also j27],|2^,|2j||3(],|in| and references 
herein. Note that the optically luminous galaxy (stellar) 
mass in clusters is about 0.19\/ft times the X-ray emit- 
ting gas mass, thus fi , = fl® n f gas (l + 0.19^/h). In what 
follows we use the f gaa values, determined by Allen et al. 
|28| from Chandra observational data, to constrain the 
parameters of the holographic dark energy model. The 
redshifts of the 26 clusters r ang e from 0.08 to 0.89. 

Following mmmmiiflim, we fit the ^ sas data 

to the holographic dark energy model described by 



(1 + 0.19^)00, 



d s A CDM (z) 



0.5 d mod 
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where d™ od and c?^ CDM are the angular diameter dis- 
tances to the clusters in the current holographic model 
and reference SCDM cosmology, respectively, and b is a 
bias factor motivated by gasdynamical simulations which 
suggest that the baryon fraction in clusters is slightly 
lower than for the Universe as a whole(see [27], and 
references herein for detailed discussions). The angular 
diameter distances to the clusters are defined as 



d A =H ( J 1 (l + z)- 1 f 
Jo 



dz' 



(11) 



where /™° d (zi;P) is computed by the holographic dark 
energy model using i|10|) . and ,fp- a ,s.i and Cf sas ,i are the 
measured value and error from |28j | for a cluster at red- 
shift Zi, respectively. The computation of \ 2 is car ~ 
ried out in a five-dimensional space, for the five param- 
eters P — (Q^ n ,c,h,fl b h 2 ,b). The probability distribu- 
tion function (likelihood) of fl® n and c is determined by 
marginalizing over the "nuisance" parameters 



c(n° m ,c) = J dhd(n° b 



h 2 )db e 



-X 2 /2 



(14) 



where the integral is over a large enough range of ft, £l°ft 2 , 
and b to include almost all the probability. We now com- 
pute £(0^,, c) on a two-dimensional grid spanned by Cl^ 
and c. The 68.3%, 95.4%, and 99.7% (namely 1, 2, and 3 
a) confidence contours consist of points where the likeli- 
hood equals e~ 2 - 31 / 2 , e - 618 / 2 j and e - 1:L - 83 / 2 f the max- 
imum value of the likelihood, respectively. 




where H 1 (here we use the natural unit, namely the 
speed of light is defined to be 1) represents the Hubble 
distance with value Hq 1 = 2997. 9ft. -1 Mpc, and E{z) = 
H(z)/H can be obtained from Q), expressed as 



E(z) 



^(l + ^) 3 xl ' 

1 - n 



X 



(12) 



Note that for the holographic dark energy model the dy- 
namical behavior of fix is determined by ©; while for 
the SCDM model we have O x = and fi° m = 1. It 
should be pointed out that the / gas data used here are 
determined assuming an SCDM model with h = 0.5. 
Hence there appears an ft/0. 5 factor in (fT77|) . We use 
the same Gaussian priors in our computation as [2^, |3(j 
with h = 0.72 ± 0.08, n° b h 2 = 0.0214 ± 0.002, and 
b = 0.824 ± 0.089, all la errors. 

To constrain the parameters of the holographic dark 
energy model, we use a x 2 statistic 
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FIG. 1: Confidence level contours of 68.3%, 95.4% and 99.7% 
in the (c, (1^,) plane. The 1 cr fit values for the parameters are: 
Q? m = 0.24tg'o5 and c = 0.6ll°;^, and the minimum value of 
\ 2 in the five-dimensional parameter space is: Xmin = 25.00. 



Figure 1 shows our main results. We plot 68.3%, 
95.4%, and 99.7% confidence level contours in the (c, $!"„) 
plane. The best fit happens at c = 0.61, f2°„ = 0.24, 
ft = 0.73, fl° b h 2 = 0.0212, and b = 0.812, with x min = 
25.00. These results are in accordance with those ob- 
tained in |3fl| where some common results, 0° = 0.24 



and Xmin ~ 25, were got from / gas fits to three models 
- ACDM model, XCDM parametrization, and cACDM 
model (quintessence with power law potential). From 
Figure 1 we see clearly that the quality of the / gas con- 
straints is much better than that of the SNe la constraints 
(see Figure 2 of 0]), namely the contours are tighter 
than those derived from SNe la data. We find, how- 
ever, that the / gas constraints on the holographic model 
are consistent with those from a joint analysis of SNe 
la, CMB, and LSS data, but the constraints from the 
latter are tighter; see Figure 6 of ^(| for comparison. 
The 1 cr fit values for c and are: c = 0.6lto'2i an( l 
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ft°„ = 0.24t°;°g. We notice that the fit value of c is less 
than 1 in 1 a range, though it can be slightly larger than 
1. This implies that according to the / gas constraints the 
holographic dark energy basically behaves as a quintom- 
type dark energy in 1 a range. 
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FIG. 2: Equation of state of dark energy w and deceleration 
parameter q, versus red-shift z, from the best fit of the / gas 
test. 
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figure, we see that the equation of state of dark energy 
w has a value of wq — —1.29 and the deceleration pa- 
rameter q has a value of go = —0.97 at present. The 
typical characteristic of the quintom-type dark energy is 
that the equation of state can cross —1. For this case, 
the crossing behavior (w(zc) = —1) occurs at a redshift 
of zc — 0.62. In addition, the transition from deceler- 
ation to acceleration (q(zx) = 0) occurs at the redshift 
zt = 0.70. Comparing our plots in Figure 2 with the 
model-independent plots in [lS|. we find that the holo- 
graphic plots for the c — 0.61 case are in good agree- 
ment with those model-independent plots for the red- 
shift range z = — 2. On the whole, the results derived 
from the / gas constraints are consistent with those from 
other cosmological tests. The parameter c which plays 
an important role in the holographic dark energy model 
is demonstrated to be less than 1 basically in 1 er range, 
which shows that the holographic dark energy tends to 
behave as quintom-type dark energy in the cosmologi- 
cal evolution. For comparing the probability distribution 
of the parameter c determined by different cosmological 
tests, we show in Figure 3 the likelihood plots of c corre- 
sponding to constraints from SNe, SNe+CMB+LSS ( for 
detail see 0), and / gas data, respectively, by further- 
more marginalizing over the "nuisance" parameter f2^. 
We see that the X-ray data provide a fairly good way for 
constraining the holographic dark energy. 

In summary, we used in this Letter the recent X-ray 
cluster gas mass fraction data from the Chandra X-Ray 
Observatory to constrain the parameters of the holo- 
graphic dark energy model. We considered a spatially 
flat FRW Universe with matter and holographic dark en- 
ergy. For the holographic dark energy model, the numer- 
ical parameter c plays a very important role in determin- 
ing the evolutionary behavior of the space-time as well as 
the ultimate fate of the Universe. The constraints from 
the /g as data show that in 1 a range the parameter c is 
basically less than 1, which implies that the holographic 
dark energy tends to behave as a quintom-type dark en- 
ergy. These constraints are consistent with those derived 
from other cosmological tests. The / gas data are proven 
to be efficacious in constraining dark energy. We hope 
that the future / gas data should provide an even tighter 
constraint on holographic dark energy model and other 
dark energy models. 



FIG. 3: Likelihood distributions of parameter c in the fits of 
SNe only, SNe+CMB+LSS, and / gas data. 



We now discuss about the cosmological consequences 
led by the best fit results of the / gas data analysis. The 
evolutions of the equation of state of dark energy and 
the deceleration parameter of the Universe correspond- 
ing to the best fit are shown in Figure 2. From this 
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